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ABSTRACT 


The stability of an infinitely long, inviscid, incompressible fluid cylinder of finite electrical 
conductivity in the presence of an axial current, a longitudinal magnetic field and an infi- 
nitely long coaxial conducting sheath placed at some distance from the fluid cylinder is in- 
_vestigated. The axial current is assumed to be maintained constant by external sources. The 
effect of surface tension is included. In the case when the fluid is mercury and the axial cur- 
rent and the longitudinal magnetic field are of the order of a few hundred amperes and a few 
hundred gauss respectively, it is shown that there exists a critical product of axial current and 
longitudinal magnetic field which determines whether the m=0 mode or the m=1 mode will 
be manifested first in the instabilities. 
The results of this theory are compared with those obtained by previous authors for the 
same geometry but with the assumption of infinite electrical conductivity for the fluid. 


I. Introduction 


_ The stability of an infinitely long, inviscid and incompressible fluid cylinder 
under the action of surface tension has been studied by Rayleigh [1]. The hydro- 
magnetic analogue of this problem (excluding surface forces) was considered, 
among other problems, by Chandrasekhar and Fermi [2] who showed that the 
presence of an.axial magnetic field frozen inside the medium increases the sta- 
bility of the cylinder for an axi-symmetric perturbation. Simon [3] generalised 
the results of Chandrasekhar and Fermi for three dimensional perturbation when 
the magnetic field is also present outside the infinitely conducting fluid cylinder. 
The extension of the stability considerations to the situation when the fluid is 
compressible and carries an axial current is of particular interest because of its 
application to the study of the behaviour of constricted gas discharges _ which 
are pinched under the mutual interaction of discharge current filaments. Kruskal 
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and Schwarzschild [4] and Taylor [5] showed that a highly constricted gas di ; 
charge is unstable. Further Taylor [6], Shafranov [7], Rosenbluth [8] and Krus 
kal and Tuck [9] showed that the stability of a pinched discharge is increaset 
(1) if it contains an axial magnetic field frozen in the medium and (2) if it 
surrounded by a coaxial conducting sheath of infinite electrical conductivity. 
In the study of hydromagnetic stability of a given configuration, it is usually 
assumed that the fluid in question has infinitely large electrical conductivity. 
This assumption is valid (1) in a highly ionised plasma and. (2) under astro- 
physical conditions. In the case of experiments carried out in the laboratory 
with mercury as the conducting fluid the hydromagnetic coupling is too weak 
to permit the assumption of infinite electrical conductivity. However, if the 
electric current is maintained by external sources at a constant value an alter- 
native simplification is allowed. On can neglect quantities proportional to the 
electrical conductivity in the induction equation and supplement the problem 
by an additional constraint on the system by demanding, for example, that the 
total current is constant when the system is perturbed. Without going into the 
question of uniqueness of the solutions obtained by such a procedure that changes 
the structure of the differential equations, it was shown recently by Murty [10] 
that the method enables one to explain satisfactorily the axi-symmetric insta- 
bilities observed by Dattner et al. [11] and Dattner [12] in their experiment which 
‘was a liquid conductor model of pinched discharges. 

In the above-mentioned experiment mercury was allowed to fall freely under 
gravity from a circular hole of 0.2 cm radius in a tank, in the form of a jet 
of 14cm length. When an axial current passes through the jet, sausage in- 
stabilities (m=O mode) are observed, and when a longitudinal magnetic field 
was also present helical instabilities (m=1 mode) were observed. 

It is the aim of the present paper to generalise the theory developed in [10] 
for a non-axisymmetric perturbation in the presence of a longitudinal magnetic 
field and also an external coaxial conducting cylinder and compare the results 
with the above-mentioned experiment of Dattner ef al. A comparison will also 
be made between the results of the present paper obtained with the assump- 
tion of vanishing electrical conductivity for the fluid with those obtained for 
infinitely conducting fluids by previous authors. 


II. Formulation of the problem 


(a) Basic assumptions and equations 


Let us consider a static, infinitely long, inviscid and incompressible fluid cyl- 
inder of radius R, immersed in a uniform magnetic field Hy, extending throughout 
space and directed parallel to the axis of the fluid cylinder. A current of con- 
stant volume density is also passing along the cylinder. The directions of cur- 
rent and the magnetic field from external sources are chosen to be parallel so 
that the resultant magnetic lines of force form right helices everywhere. Let 
the fluid be non-magnetic, have vanishingly small electrical conductivity and be 
surrounded by an infinitely conducting coaxial cylinder placed at a radial dis- 
tance equal to AR, from the axis. We shall further assume that the region 
between the fluid cylinder and coaxial conducting cylinder is a vacuum. The 
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I= constant. 


ne symbols in eqns. (2.1) to (2.6) have their usual meaning [13]. Eqn. (2.7) 
is a constraint on the system in the present case which indicates that the axial — 
current in the fluid cylinder is maintained constant by external sources, 


« _ 


‘ (b) Equilibrium 


The equilibrium state of the fluid cylinder is described by 


a 0=—V pet+jex H., 

and eqns. (2.3) and (2.7). If we choose a system of right circular cylindrical 

coordinates with the z-axis coinciding with the axis of the fluid cylinder and 
the current density inside the fluid is given by 


= . . 

j je=(0, 0, Jo), (2.8) 
4 - 

the fluid cylinder will be in equilibrium with the pressure distribution given by 
1 The inclusion of the gravity can have a nontrivial effect on the solution of this stability 
_ problem only when fluid velocity is a function of the position along the axis of the fluid cylinder. 
a In the case of Dattner’s experiment there is in fact a variation of the fluid velocity along the 
axis of the fluid cylinder and the equilibrium distribution of the physical quantities like the 
current density, the pressure and the magnetic field are functions of not only the 7-coordinate 
- but also of the z-coordinate. Since the study of stability of such an equilibrium presents enormous 
difficulties, we shall study only an idealised model by assuming that the fluid cylinder is in static 
equilibrium defined by functions depending only on the y-coordinate. We shall see that even 
- this simplified model of Dattner’s experiment is not so bad as is might appear at first sight. 
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De=nje(Ri—r)+T/Ry, when r<Rp, arb’ 
where 7’= surface tension of the fluid. The pressure outside the fluid is ca 


because of the basic assumptions we have made. | 


(c) Perturbation and criterion of instability 


In order to determine whether the equilibrium state given by eqns. (2.8) and | 
(2.9) is stable or not, we shall assume that the fluid cylinder is subjected to 
a perturbation so that its deformed surface is given in circular cylindrical co- 
ordinates by 


r=R+a cos (—m@0+kz+<8), (2.10) 


where ¢ is a phase factor and a may be a function subject to the condition 
|a|/R<1, so that we can neglect quantities ~(a/R)” in the linear approxima- 
tion. The basis for the assumption regarding the form of the perturbed sur- 
face as given by eqn. (2.10) is that a general form of the perturbed boundary 
of the fluid cylinder can be obtained by a combination of such terms as long 
as nonlinear effects are ignored. It is therefore sufficient to discuss the effect 
of one typical term as the one given in eqn. (2.10). It is no loss of generality 
in the present discussion of stability of the equilibrium to assume that the phase 
factor e« in eqn. (2.10) is equal to zero. We further assume that a is periodic 
in time so that 


r=Rh-+a‘cos 9, (2.11) 
and a=, cos wt, (2.12) 
where @o= — m0. kz. (2.11 a) 


The purpose of the stability theory is to calculate the values of wm. If w is 
found to have an imaginary value for any value of m and k, we conclude that 
the fluid cylinder is unstable for arbitrarily small perturbations of the boundary. 

Since the fluid is assumed to be incompressible one can show that R and R, 
are related by 


Ro= R? + a?/2. (2.13) 


Eqn. (2.13) shows that in the linear approximation one can ignore the differ- 
ence between F and Rp. 


III. Solution for @ 


(a) Current density distribution 


If we combine eqns, (2.4) and (2.5) we obtain 
oH ue f Wl 
—= —Vx|*-—vx 
Kom > (3.1) 
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all other quantities in eqn. (3-1) are finite. The current density vector 
ilso solenoidal (see eqn. (2.3)) and therefore eqn. (3.2) shows that we can 
a scalar function y such that 


; j=2f"t+avy, 148.8) 
where y satisfies v* p=, ce 
d j* is the mean current density flowing in the z-direction. The solutions of 


n. (3.4) (in circular cylindrical coordinates) that are regular in the domain 


secupied by the fluid can be written as . 
y=CI,, (kr) sin 9, (3.5). 


where J, (kr) is the modified Bessel function and C is an arbitrary constant. We note 
e that the eqn. (3.3) is formally the same as the one obtained by Chandra- 
ekhar and Fermi [2] for the distribution of the magnetic field in connection 
th the study of the stability of a fluid cylinder in the presence of an axial 
‘magnetic field frozen in the medium, and the expression for wy in the eqn. (3.5) 
is a generalisation of their eqn. (75) in [2]. Since we are interested only in the 
“Tinear approximation, we have neglected higher order terms in eqn. (3.5). 

_ The boundary condition for the solution of eqn. (3.4) is that the normal 
-eomponent of the current density vector at the perturbed boundary is zero, or? 


(j-D)o=9, (3.6) 


where n is the unit vector normal to the perturbed boundary and the suffix b 
indicates that the quantity concerned shall be evaluated at the perturbed bound- 
ary i.e. at r=R+acosqg. One can show that n is given by 


n=(l1, 0, 0)+(0, —am/R, ak) sin 9. (3.7) 


“Using eqns. (3.3), (3.5), (3.6) and (3.7) we find that 
C= —j*/In (ER), (3.8) 
where (and in the following) the prime to the modified Bessel function indicates 


differentiation with respect to the argument. Due to the constraint that the 
total current flowing along the fluid cylinder is constant (eqn. (2.7)), we find that 


J” =o. (3.9) 


1 The eqn. (4.6) of [10] should read as (j-M)r=R+a cos ke = 0. 
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We note that if we put m=0 in eqns. (3.11) we obtain eqns. (4.8) of 

: miljurdmaahl3 104) -«Ageroere eagerly ” 
(b) Magnetic field distribution = Lgdeg 


In order to obtain the perturbed magnetic field distribution, H,, inside t 24 
fluid cylinder from eqns. (2:3) and (2.6) we assume that we can separate 
variables r, @ and z, and write 


HH; = (A, ;, Agi, Ia), 


2% 


. 


Pa 
it 


(3.12) 


at 


where H,; = H,, (r) sin q, x . 
Ha = He (r) cos p+ 2a jor, ~ (3.13) 
and . Hy=H,,(r) cos ¢. 


The phases in eqns. (3.13) for the magnetic field are so chosen as to be con- 
sistent with the phases of the current density distribution given by eqns. (3.11). 
The reason for the particular choice of the 6-component of the magnetic field 
in eqns. (3.13) is that the z-component of the current density has a constant 
mean value which gives, inside the fluid cylinder, a toroidal component of the 
magnetic field increasing proportionally with the distance from the axis. If we 
eliminate the r- and §-components of the magnetic field between the eqns. (2.3) 
and (2.6) we find, using eqns. (3.11), that H,, (r) satisfies the equation 


|? < (: ra — ( + "| H.(r)=0, (3.14) 


whose well-behaved solution in the domain of interest can be written as 
H, (r) =A, (kr), (3.15) 

where A is an arbitrary constant. Using eqns. (2.3), (3.11), (3.13) and (3.15), 

one can show that the radial parts of the remaining components of the magnetic 


field inside the fluid cylinder are given by 
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Brg: 474M jy Im (kr) 
BASAL. (bie 
(7) ( 8 5 Hs gelae 
(3.16) 
mAI,,(kr) 42a 
H 24 sige bell er pe Jo 7" 
6 (7) BF I, (eR) J pw AkY) 


z 

~In the region outside the fluid cylinder and inside the coaxial conducting 
cylinder the magnetic field vector is both irrotational and solenoidal and hence 
can be obtained as the gradient of a harmonic function. One can show that 
the; general form of the magnetic field distribution, Hp, outside the fluid cyl- 
inder is given by 


Hy =[H,9, Hoo, Heol; (3.17) 
where A, 9=k[B, In (kr) + By Kn (kr)] sin 9, 
m 224, R? 
Hog = ~~ [By In (kr) + By Km (Kr)] c08 oo oo (3.18) 
and A. ,=k (By Im (kr) + By Kn (kr)] cos 9, 


where) K,, (x) is that modified Bessel function that is regular at 4 = 0 and B, 
and B, are arbitrary constants. The magnetic field should be represented by 
real functions and therefore we assume that k>0 in eqns. (3.18). Since the 
external cylinder at a radial distance equal to A R, is assumed to be an ideal 
conductor, the radial component, H,,, of the magnetic field must vanish at 
r=A R,, which shows that 

B, Kn (AkR,) 


B, 74 ie R,) ad) 

We assumed that the fluid cylinder has only a volume distribution of cur- 
rent and therefore the magnetic field vector has to be continuous at its boundary. 
By matching at the boundary the inside and outside solutions of the magnetic 
field as given by eqns. (3.13) to (3.19), we find that the constant A in eqns. 
(3.15) and (3.16) is given by 


470M jo Ln (kB) 


TER) (3.20) 


Ax — 


K,(AkR 
Kn (i R) — In ame ee, 


where the difference between R and R, has been ignored by virtue of eqn. (2.13). 
The constants B, and B, can similarly be determined but their values are not 
given since they are not required for our purpose. Therefore eqns. (3.13) to 
(3.20) determine the perturbed magnetic field distribution both inside and out- 
side the fluid cylinder up to terms of order a. 


(c) Pressure distribution 


If we take the divergence of both sides of eqn. (2.1) und use eqns. (2.2), 
(2.3), (3.3) and (3.11) we obtain, ignoring higher order terms, 
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Ly (er) 3.21 
I, (eR) °°” oa 


V2 p= —4nje+8nak jo 


t 
which is merely Poisson’s equation. Following the method adapted in [10] one 
can show that the solution of eqn. (3.21) is given by 


fl “ A r Im (kr) 
Py ies) Vee ay cos p 
sgl hg le OE 2 “2 | 

He - |= (ke —1)-2 R i 3.229 | 
a7) | AR +m 1)—22 70 R| cos p ( | 


Eqn. (3.22) represents the pressure distribution at all points in the interior of 
the perturbed boundary of the fluid cylinder. It has been obtained from the 
condition that the pressure at a point just inside the boundary is that due to 
surface tension. 


(d) Dispersion relation 


In order to obtain the values of w, we note that the radial part of eqn. (2.1) 
evaluated at the boundary can be written, using eqns. (2.11) and (2.12), as 


—o oa cos p={—8 p/ér+ (jxB),}», (3.23) 


where the suffix to the braces indicates that the quantity concerned must be 
evaluated at r=R+a cosy. Using eqns. (3.11), (3.13), (3.15) to (3.20) and (3.22) 
and including in the electromagnetic force the effect of the homogeneous longi- 
tudinal magnetic field produced by external sources, one can show that eqn. 
(3.23) leads to 


T x). 47? m I Hy Im (#) 
CN 2 ra i eee aes 0 x 
o R® fs (a) (a +m ) zo R* m (x) zo R® i. (ar)’ 


where [=7j, R®, x=kR and 


8 
at 
= 
& 
“I 
~) 


@ 


(3.24) 


_, , %fIn(e) In (e))  m? Ip (x) od is Im (%) Km (A 2) 
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(3.25) 


and H,—= homogeneous magnetic field parallel to the z-axis produced by external 
sources. Hqn. (3.24) is the required dispersion relation which gives the angular 
“frequency” of vibration of the surface of the fluid cylinder for various modes 
of perturbation given by eqn. (2.11). 


IV. Discussion 


(a) General results 


As was stated in § II (c), the fluid cylinder is unstable for arbitrarily small 
perturbations for all values of x for which w? assumes a negative value. In 
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_ present solution given by eqn. (3.24) only real values of x>0 are allowed 
nd w” is then a real number showing that there are no states of over-stability, 
result which is similar to the one obtained by Bernstein et al. [14]. If we set 
)=0 in eqn. (3.24), we obtain the eqn. (6.2) of [10]. When m=0, the effect 
f an external magnefic field, Hy, is zero, which is a consequence of the fact 
that we have assumed the fluid to be a poor conductor. When m+0, we find 
. only the third term on the right-hand side of eqn. (3.24), containing the 
contribution due to an external magnetic field, changes sign if m or a changes 
= In the present problem we assumed that the magnetic field produced 
by external sources is directed parallel to the direction of axial current, and 
so the fluid cylinder is unstable for that mode m0 which is right-handed in 
a right circular cylindrical coordinate system. An examination of the current 
distribution given by eqns. (3.11) shows that the unstable deformation is the 
one that makes the current density distribution tend to the force-free state 
‘at the points that are perturbed away from the axis. One can also verify that 
such a distribution of current increases the magnetic pressure at a point outside 
that part of the fluid surface which is perturbed inwards and vice versa and 
one concludes that the cylinder is unstable for such perturbations. We note here 
that the first term on the right-hand side of eqn. (3.24) is the one obtained by 
Rayleigh [1]. 


t 


(b) Neither external conducting cylinder nor surface tension 


If the jet is not enclosed in a conducting cylinder, i.e. if A= oo and the fluid 
‘has no surface tension, one can show that eqn. (3.24) reduces to 


és 42 [© {Iin(z) Im (x)| m? Ip (2) eee I, (x) 
@ = wo Rt E 2 ir (x) Ii, ae ed, (#) {1+ Im (20) Km (#)} wok Ii, (x) 
(4.1) 


where the suffix 6 in wz indicates the conditions stated in the beginning of this 
subsection for which the equation is valid. By considering the asymptotic be- 
haviour of the modified Bessel functions one can show 


2 
Lit wi> ms z (|m|—1)+ 0 (2), (4.2) 
o+0 wok 
ap be m RH ‘) 
2 ps , 4.3 
and peli > ale Wi o(: (4.3) 


> 


The behaviour of w? for «0 shows that for infinitely long wavelength dis- 
turbances the fluid cylinder is neutral for the |m|=1 mode, while it is stable 
for |m|>1. In this respect the effect of an axial current is similar to that of 
surface tension [1]. The behaviour of w2 for x>co shows that the fluid cyl- 
inder is unstable for a disturbance of very small wavelength, but the growth 
rate of perturbations is finite for all finite values of m, which is in contrast 
with the infinite growth rate of small wavelength disturbances obtained by Krus- 
kal and Schwarzschild [4] and Taylor [5], who assumed surface current distri- 
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bution for an infinitely conducting plasma cylinder. Also, : Taylor obt in 
finite growth rate for small wavelength disturbances if the infinitely conduc 
plasma has a volume current distribution [5]. We see therefore that the growth 
rate is finite for volume distribution of current and is infinite for surface cur- 
rent distribution. The reason for this behaviour is that the effective mass in- 
volved in the motion decreases with wavelength of perturbation, while the sur- 
face density of the force remains constant so that the acceleration becomes very 
large in the limit of vanishingly small wavelengths. It must, however, be noted 
that the fluid equations and the perturbation technique are no longer valid | 
the limit of vanishing wavelength of the perturbation. : 

As mentioned earlier, the asymptotic values of w;, given by eqns. (4.2) and 
(4.3), show that all the modes are unstable for small wavelengths, while for 
long wavelengths the modes with |m|>1 are stable. Therefore a given mode 
|m|>1 passes through a neutral state as the wavelength of the perturbation 
decreases from infinitely large values to small values. The variation of w5 for 
m=2 when H,=0 is shown in Fig. 1, where it is seen that «> has a zero when 
x4, The neutral point for a general value of m>1 can be obtained by con- 
sidering the asymptotic behaviour of the modified Bessel functions in the range 
1<uz<co, One can show that in the above limits and when H,=0, w5=0 
when mw 2/2. 

The reason for the existence of a neutral state for some wavelength when 
m=2 is qualitatively as follows. An examination of the current density distri- 
bution given by eqns. (3.11) shows that for deformations of finite wavelength 
and m=2 there are both parallel components (z-components) and anti-parallel 
components (#-components) of the current. The former components attract each 
other, while the latter repel, so that there can be a combination of these cur- 
rent configurations with the net force equal to zero. The situation for m=2 is 
shown schematically in Fig. 2. 

It is interesting to compare the above results with those obtained by Taylor 
[5], who showed that a compressible, infinitely conducting fluid cylinder with 
surface current distribution is neutral to a perturbation of a wavelength corre- 
sponding to #~m* when m>1. He has also shown that such a neutral state 
does not exist when the fluid has a volume current distribution. In fact the 
fluid cylinder is completely stable for all the modes with m>1 when the cur- 
rent is evenly distributed in the whole cross-section of the cylinder and this 
result is independent of the compressibility of the medium. We see therefore 
that the behaviour of the modes with m>1 is strongly dependent on the con- 
ductivity of the fluid. Although the volume distribution of the current gives 
greater stability to the (infinitely conducting) fluid cylinder than a surface distri- 
bution of current (in the absence of frozen axial magnetic field), the modes with 
m>1 are unstable if the conductivity of the fluid is small and the currents are 
maintained by external sources while they are completely stable if the con- 
ductivity is infinite as considered by Taylor. 

An examination of the limiting form of w? when 20 (eqn. 4.2)) shows that 
when m=2 the frequency of vibration, y=«,/22, of the surface is given by 


PAD epee es 
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Fig. Fe Variation of 7@ R4a5/I? for the mode m= 2 with «=k R. The points above the abscissa 
indicate oscillation of the fluid cylinder for the corresponding values of « with the mode m= 2. 
The point: where the curve cuts the abscissa corresponds to a critical wavelength for which the 
ae fluid cylinder is neutral to a perturbation of that wavelength of the mode m= 2. 
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Fig. 2. Effect of a perturbation of finite wavelength with m= 2 without any external magnetic 
field. The figure shows qualitatively the cross-section of the perturbed fluid cylinder. The cur- 
rent distribution has parallel components (z-components), whose mutual attraction is denoted 
by F, and eh the anti-parallel components (0-components) repel each other, whose effect is de- 


noted by #, and 2a These forces can balance each other for a certain wavelength of the per- 


turbation so that the fluid cylinder is neutral for that wavelength of the perturbation. 
i 


where H,(R) is the toroidal component of the magnetic field at the boundary 
of the cylinder due to the axial current. Eqn. (4.4) suggests that the frequency, 
y, is roughly equal to the inverse of the time of travel of an Alfvén wave 
around the circumference of the fluid cylinder. It is also interesting to note from 
eqn. (4.3) that, when zoo, ie. for a disturbance of very small wavelength, 
the e-folding time of the disturbance is also of the same order as the time of travel 
of an Alfvén wave around the circumference a result which is similar to the 
one obtained by Taylor [5] for an infinitely conducting fluid, the only difference 
being that in his case the e-folding time of the disturbance is /2 times the one 
obtained by us in the above problem. We may conclude from the above results 
that the time scales defined by Alfvén wave propatation over characteristic dis- 
tances of the problem play a significant role in both infinitely conducting fluids 
as well as finitely conducting fluids, provided the currents are maintained by 
external sources in the latter case. 
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(c) No external conducting cylinder 


Tf we take into account the effect of surface tension of the fluid but assume 
that A= oo, we get from eqn. (3.24), writing w= We 


T x2 Lin (x) j ' 2 
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Fig. 3. Fig. 4. 


Fig. 3. Variation of x, with H,, the homogeneous magnetic field produced by external sources 
nea two values of the axial current flowing in a mercury cylinder of 0.2 cm radius. The quanti ; 
2, is proportional to the wave number for which the fluid cylinder is neutral to the perturba- 
tion of the mode m=1. The neutral state arises because of the stabilising action of surface ten- 
sion of the fluid. 


Fig. 4. Variation of x,, with H,, the homogeneous magnetic field produced by external sources 
for two values of axial current flowing in a mercury cylinder of 0.2 cm radius. The quantity 
. &, is proportional to the wave number of the perturbation of the m=1 mode for which the fluic 
cylinder has a maximum growth rate. The existence of the wavenumber with a maximum growt 
rate for perturbation of the mode m=1 is due to the presence of surface tension of the fluid 
(mercury in the present case). 


where w; is given by eqn. (4.1). The first term in eqn. (4.5) is due to the 
_ effect of surface tension which has a stabilising effect for all wavelengths when 
m>1 and can give rise to instability for some range of wavelengths (less than 
the circumference of the fluid cylinder) when m= 0, a result which was discussed 
by Rayleigh [1]. When zoo, that is, when the wavelength tends to zero, the 
term in eqn. (4.5) containing the effect of surface tension tends to infinity for 
all values of m. The stabilising effect of surface tension at small wavelengths 
is similar to that of an axial magnetic field frozen into an infinitly conducting 
fluid cylinder [6, 7, 8, and 9]. When m=1 and a->0, one can show that w2—+>0 
and is negative for small values of x. Since w2>0 for a—>co, due to the effect 
of surface tension, we see that w2 has a zero for some value of 2 which we 
may call 2 and has a minimum at some value of x which we may call ap. 
(The suffix m in x, should not be confused with the same symbol denoting the 
mode of perturbation.) The values of 2, and 2,, for various values of axial 
current and longitudinal magnetic field for the mercury cylinder of 0.2 cm radius 
for the mode m=1 are obtained by a graphical method and are shown in 
Figs. 3 and 4. We see from these figures that both «, and z,, increase with 
the current and the magnetic field which means the limit of stability moves 
towards smaller wavelengths with increasing values of current and magnetic field. 
The ratio «,/«,, is approximately constant equal to 1.9 in the range of the values 
of the current and the magnetic field shown in Figs. 3 and 4. 

Due to the stabilising effect of surface tension all modes with m>1 are stable 
for the values of the currents and the magnetic fields used in the experiment 
of Dattner [12] and the only modes that are not stabilised are those with m=0 
and 1. The behaviour of the mode m=0 has been discussed in [10]. For m=] 
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. Variation of the maximum growth rate, | ,,|, (in non-dimensional units) for the m=1 
of perturbation of the mercury cylinder ot 0.2 cm radius with the homogeneous magnetic 
produced by external sources for two values of the axial current. @ is the density of the 
4 fluid (mercury); 7’ is the surface tension; R is the radius of the cylinder. 


non-dimensional maximum growth rate,! (9 R°/T)?|@,|, corresponding to a, of 
. 4 for given values of axial current and longitudinal magnetic field for the 
reury cylinder of 0.2 cm radius (and 7’=487 dyne cm~’ and 9 =13.6 g cm~*) 
shown in Fig. 5. The experimental value of the same nondimensional growth 
for the m=1 mode which is seen to develop in the experiment of Dattner 
11, 12] with a longitudinal magnetic field of 300 gauss and an axial current 
of 100 A is 1.8+0.6 whose mean value is about twice larger than the theo- 
tetical value shown in Fig. 5. The experimental value is obtained from the 
BRcsromonts on deformations which can by no means be considered as small 
‘perturbations of the fluid cylinder. Taking into consideration this fact and the 
possibile errors in actual value of the surface tension which can be smaller than 
that employed in the theoretical calculation the agreement between theory and experi- 
ment within a factor of 2 is not unreasonable. The observed large value of the 
growth rate may be due to non-linear effects but since there are other sources 
‘of error due to gravity and viscosity which are not included in the theory it 
is difficult at this stage to give a reasonable explanation for the discrepancy 
between theoretical and experimental value. 
~ 
; (d) Relative growth rates of the m-0 and the m=1 modes 
- We have seen before that the m=0 mode is unaffected by the presence of 
a longitudinal magnetic field, while modes m1 can be made unstable. When 
the effects of the surface tension and the longitudinal magnetic field are in- 
cluded, the growth rate for the m=1 mode can be smaller than that for the 
m=0 mode for some range of values of the axial current and the longitudinal 
magnetic field due to the stabilising action of the surface tension for the m=1 
mode. If the surface tension were zero, only the modes with m>1 will develop 
for all non-zero values of the longitudinal magnetic field. In Fig. 6 are shown 
the growth rates for the m—=0 mode and the m=1 mode for various values of 
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c. s. muRTY, Instability of a conducting fluid cylinder 
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Fig. 6. Fig. 7. 


Fig. 6. Relative growth rates of the m=0 mode and the m=1 mode of perturbation of the 
mercury cylinder of radius 0.2 cm. The ordinate represents the growth rate in non-dimensional 
units and the abscissa represents the axial current in the cylinder. The continuous curves re- 
present the growth rate of the m=1 mode for various values of the homogeneous magnetic field 
produced by external sources and the discontinuous curve represents the same quantity for the 
m=Q0 of mode of perturbation, (Note that the m=0 mode is unaffected by the homogeneous 
magnetic field.) At the intersection of the continuous curves with the discontinuous curve the 
values of the growth rates for the m= 1 mode and the m=0 mode are equal, The corresponding 
values of the axial current and homogeneous magnetic field make a product approximately 
equal to 10* A gauss. In the range of the values of the current and magnetic field shown in the 
figure the mercury cylinder is stable for all the modes with m> 1 due to the stabilising action of 
the surface tension. 


Fig. 7. Effect of an infinitely conducting coaxial cylindrical conductor enclosing the mercury 
cylinder of 0.2 cm radius. The ordinate represents the square of the “frequency” (in non-dimen- 
sional units) of “‘oscillation” of the fluid cylinder for perturbation of the mode m=1 and the 
abscissa represents the corresponding wavenumber (in non-dimensional units) of the perturba- 
tion. The axial current is 100 A and the homogeneous magnetic field is 300 gauss. The effect 
of surface tension of mercury is included. The values of the parameter A are equal to the ratio 
of the radial distance of the infinitely conducting coaxial cylinder to the radius of the mercury 
eylinder. The points below the abscissa correspond to instability of the mercury cylinder at the 
respective wave numbers for the m=1 mode of the perturbation. The stabilising effect of the 
conducting sheath decreases rapidly with increasing values of A. 


the current and the longitudinal magnetic field. We see from this figure that 
for a given value of the longitudinal magnetic field, say 100 gauss, the growth 
rate for the m=1 mode is smaller than that for the m=0 mode up to the value 
of the axial current of approximately 100 A. It is seen from the curves in 
Fig. 6 that the growth rates of the m=0 mode and the m=1 mode are equal 
when the product of the axial current and the longitudinal magnetic field is 
approximately equal to 10* A gauss. In the experiment of Dattner [12] the values 
of the current and the longitudinal magnetic field used give a product which 
is larger than the above value. It is of some interest to repeat the experiment 
to verify the “‘critical product’? of the axial current and the longitudinal mag- 
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(e) Effect of an external conducting winder 


f the fluid cylinder is surrounded by an infinitely conducting coaxial cylinder, 
stability of the fluid cylinder is governed by eqn. (3.24). The value of 
/T)o" is plotted versus « in Fig. 7 for m=1, taking H,—300 gauss and 
00 A for a mercury cylinder of 0.2 cm radius for three values of A, namely, 
2.0 and co. (Note that the m+0 modes only are affected by the presence 
mn external conductor.) We see from this figure that the external coaxial 
uctor has a stabilising effect but cannot completely stabilise the fluid cyl- 

for the above values of the axial current and the longitudinal magnetic 
Id even when placed at a radial distance of 1.2 times the radius of the fluid 
ylinder. The stabilising effect decreases very rapidly with increasing radial dis- 
ce of the conductor. The result should be compared with the corresponding 
ization in an infinitely conducting plasma cylinder which is shown to be stab- 


ie ano 


e absence of such a marked stabilising effect in our present example is due 
the assumption of poor electrical conductivity for the fluid. 
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